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The theoretical study 01 HI] and experimental ob- 
servation ]4J of the quantum spin Hall effect in 2D sys- 
tems, followed by the discovery of a similar phenomenon 
is 3 dimensions JUSlTlSlit], have generated consider- 
able interest in topological states of free electrons. Both 
kinds of systems are time-reversal invariant insulators. 
More specifically, they consist of (almost) noninteract- 
ing fermions with a gapped energy spectrum and have 
both the time-reversal symmetry (T) and a f/(l) symme- 
try (0. The latter is related to the particle number, which 
is conserved in insulators but not in superconductors or 
superfiuids. Topological phases with only one of those 
symmetries, or none, are also known. Such phases gen- 
erally carry some gapless modes at the boundary^ 

The classification of gapped free-fermion systems de- 
pends on the symmetry and spatial dimension. For exam- 
ple, two-dimensional insulators without T symmetry are 
characterized by an integer v, the quantized Hall conduc- 
tivity in units of e 2 /h. For systems with discrete trans- 
lational symmetry, it can be expressed in terms of the 
band structure (more exactly, the electron eigenstates as 
a function of momentum); such an expression is known 
as the TKNN invariant 111 ill , or the first Chern number. 
A similar topological invariant (the k-th Chern number) 
can be defined for any even dimension d. For d = 0, it is 
simply the number of single-particle states with negative 
energy (E < Ef = 0), which are filled with electrons. 

However, the other three symmetry types (no symme- 
try, T only, or both T and Q) do not exhibit such a simple 
pattern. Let us consider systems with no symmetry at all. 
For d = 0, there is a Z2 invariant: the number of elec- 



1 In contrast, strongly correlated topological phases (with anyons in the 
bulk) may not have gapless boundary modes [10]. 



trons (mod 2) in the ground state. For d = l,a system in 
this symmetry class, dubbed "Majorana chain", also has 
a Z2 invariant, which indicates the presence of unpaired 
Majorana modes at the ends of the chain 11211 . But for 
d = 2 (e.g., a p x + ip y superconductor), the topological 
number is an integer though an even-odd effect is also 
important fiH[l4ll . 

T-invariant insulators have an integer invariant (the 
number of particle-occupied Kramers doublet states) for 
d = 0, no invariant for d = 1, and a Z2 invariant for 
d = 2 H H and for d = 3 fl 0, 0]. 3D crystals (i.e., 
systems with discrete translational symmetry) have an 
additional 3Z2 invariant, which distinguishes so-called 
"weak topological insulators". 

With the exception just mentioned, the topological 
numbers are insensitive to disorder and can even be de- 
fined without the spectral gap assumption, provided the 
eigenstates are localized. This result has been established 
rigorously for integer quantum Hall systems [! 15lll6lll7ll . 
where the invariant v is related to the index theory and 
can be expressed as a trace of a certain infinite operator, 
which represents the insertion of a magnetic flux quan- 
tum at an arbitrary point. Its trace can be calculated with 
sufficient precision by examining an /-neighborhood of 
that point, where I is the localization length. A similar lo- 
cal expression for the Z2 invariant of a ID system with no 
symmetry has been derived in Appendix C of Ref. [ 14]; 
it involves an infinite Pfaffian or determinant. 

In this paper, we do not look for analytic formulas for 
topological numbers, but rather enumerate all possible 
phases. Two Hamiltonians belong to the same phase if 
they can be continuously transformed one to the other 
while maintaining the energy gap or localization; we will 
elaborate on that later. The identity of a phase can be de- 
termined by some local probe. In particular, the Hamil- 



TABLE 1. Classification of free-fermion phases with all possible combinations of the particle number conservation (Q) and 
time-reversal symmetry (T). The 7to(C q ) and 7to(R q ) columns indicate the range of topological invariant. Examples of topologically 
nontrivial phases are shown in parentheses. 
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tonian around a given point may be represented (in some 
non-canonical way) by a mass term that anticommutes 
with a certain Dirac operator; the problem is thus reduced 
to the classification of such mass terms. 

Prior to this work, there have been several results to- 
ward unified classification of free-fermion phases. Alt- 
land and Zirnbauer 111 811 identified 10 symmetry classes 
of matrices^ which can be used to build a free-fermion 
Hamiltonian as a second-order form in the annihilation 
and creation operators, dj and a}. The combinations of 
T and Q make 4 out of 10 possibilities. However, the 
symmetry alone is only sufficient to classify systems in 
dimension 0. For d = 1, one may consider a zero mode 
at the boundary and check whether the degeneracy is 
stable to perturbations. For example, an unpaired Majo- 
rana mode is stable. In higher dimensions, one may de- 
scribe the boundary mode by a Dirac operator and like- 
wise study its stability. This kind of analysis has been 
performed on a case-by-case basis and brought to com- 
pletion in a recent paper by Schnyder, Ryu, Furusaki, and 
Ludwig 11911 . Thus, all phases up to d = 3 have been char- 
acterized, but the collection of results appears irregular. 

A certain periodic pattern for Z 2 topological insula- 
tors has been discovered by Qi, Hughes, and Zhang H2Qfl . 
They use a Chern-Simons action in an extended space, 
which includes the space-time coordinates and some pa- 
rameters. This approach suggests some operational inter- 
pretation of topological invariants and may even work for 
interacting systems, though this possibility has not been 
explored. In addition, the authors mention Clifford alge- 
bras, which play a key role in the present paper. 



2 These classes are often associated with random matrix ensembles, but 
the symmetry pertains to concrete matrices rather than the probability 
measure. 



We report a general classification scheme for gapped 
free-fermion phases in all dimensions, see Table Q] It ac- 
tually consists of two tables. The small one means to 
represent the aforementioned alternation in TR-broken 
insulators (a unique trivial phase for odd d vs. an inte- 
ger invariant for even d). The large table shows a pe- 
riod 8 pattern for the other three combinations of T and 
Q. Note that phases with the same symmetry line up di- 
agonally, i.e., an increase in d corresponds to a step up 
(mod8). (^-invariant ID superconductors were studied 
in Ref. 12111 . The (p x +iPv)^ + (Px—ipy)l phase was pro- 
posed in Refs. ||22[ l23i fl9ll : the last paper also describes 
an integer invariant for 3 He-B.) The 2 + 8 rows (in- 
dexed by q) may be identified with the Altland-Zirnbauer 
classes arranged in a certain order; they correspond to 2 
types of complex Clifford algebras and 8 types of real 
Clifford algebras. Each type has an associated classify- 
ing space C q or R q , see Table [2] Connected components 
of that space (i.e., elements of 7Co(Rq) or 7to(C q )) corre- 
spond to different phases. But higher homotopy groups 
also have physical meaning. For example, the theory pre- 
dicts that ID defects in a 3D TR-broken insulator are 
classified by %\ (C\ ) = Z. 

The (mod2) and (mod8) patterns mentioned above 
are known as Bott periodicity; they are part of the math- 
ematical subject called K-theory. It has been applied 
in string theory but not so much in condensed matter 
physics. One exception is Hofava's work [24] on the clas- 
sification of stable gapless spectra, i.e., Fermi surfaces, 
fines, and points. In this paper, we mostly use results 
from chapters II— III of Karoubi's book 1I25I1 . in particular, 
the relation between the homotopy-theoretic and Clifford 
algebra versions of ^-groups (a variant of the Atiyah- 
Bott-Shapiro construction II26I1 ). 



TABLE 2. Bott periodicity in complex and real -theory. (The parameters k,m,n should be taken to infinity.) 
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Above: The classifying space Co parametrizes Hermitian ma- 
trices X with ±1 eigenvalues. C q is the 17-th loop space of Cq; 
it parametrizes such matrices X that anticommute with q Clif- 
ford generators. 

Right: Similar classification for real symmetric matrices. 
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SOME EXAMPLES 

To get a glimpse of the mathematical structure underly- 
ing the topological classification, we consider a second- 
order transition between two phases, where the energy 
gap vanishes at some value of parameters. In this case, 
the low-energy Fermi modes typically have a Dirac spec- 
trum, and the phases differ by the sign of the mass term. 

Let us begin with the simplest example, the Majorana 
chain |12||. This model has one spinless Fermi mode per 
site, but the number of particles is not conserved, which 
calls for the use of Majorana operators: 



C2/-1 =a r 



J' 



(; = !,...,«). (1) 



By convention, operators acting in the the Fock space (as 
opposed to the mode space) are marked with a hat. The 
Majorana operators are Hermitian and satisfy the com- 
mutation relations c/c m + c,„c/ = 25/ m ; thus, c\,... ,C2„ 
may be treated on equal footing. (But it is still good to 
remember that 621-1 and C2j belong to the same site j.) 
The advantage of the Majorana representation is that all 
model parameters are real numbers. 

A general free-fermion Hamiltonian for non- 
conserved particles has this form: 



* j,k 



(2) 



where A is a real skew-symmetric matrix of size 2n. The 
concrete model is this: 



n-\ 



M I^2/_iC2/ + v£c2/C2M-i . (3) 

z V 1=1 1=1 I 



At the transition between 'the "trivial phase" (\u\ > |v|) 
and the "topological phase" (\u\ < |v|), there are two 
counterpropagating gapless modes. They may be repre- 
sented by two continuous sets of Majorana operators, 
fjj (x) (j = 1,2). The effective Hamiltonian near the tran- 
sition point has this form: 



H 2.1 77 \-m -d 



fj dx, 



1 



(4) 



where m ~ u — v. Thus, we need to study the Dirac 
operator D = yd + M, where 7 = a z and M = mia y . 
If m gradually varies in space and changes sign, e.g., 
m(x) — —ax, the Dirac operator has a localized null 
state, which corresponds to an unpaired Majorana mode 
in the second quantization picture. The existence of the 
true null state is a subtle property, but it has a simple 
semiclassical analogue: a continuous transition between 
a positive and a negative value of m is impossible without 
closing the gap. 

We now consider a model with two real fermions prop- 
agating in each direction, so that the mass term has more 
freedom. This situation occurs, for example, at the edge 
of a 2D topological insulator. A gap opens in a magnetic 
field or in close contact with a superconductor 02711 . The 
Hamiltonian is as follows: 



/ i// T + i// T f \ 



H 



7 = 



f\ T (yd +M)fjdx, fj 
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-h z {iay) 



m 

h z (ia y ) 



(5) 



(6) 



m = -h x {io y ) + h y I - (Re A)a x - (ImA)cr\ (7) 

If h z — 0, the energy gap is given by the smallest singu- 
lar value of to; it vanishes at the transition between the 
"magnetic" and "superconducting" phase as the function 
det(ro) = h% + hZ — |A| 2 passes through zero. The pres- 
ence of h z complicates the matter, but if the spectrum 
is gapped, h z can be continuously tuned to zero without 
closing the gap. We will see that, in general, the mass 
term can be tuned to anticommute with 7, in which case 
M consists of two off-diagonal blocks, m and — m T . 

With n modes propagating in each direction, the non- 
degenerate anticommuting mass term is given by m 6 
GL(n,R). This set has two connected components, hence 
there are two distinct phases. Note that the set GL(«,R) 
is homotopy equivalent to R\ = 0(n) (see Table |2); it 
provides the classification of systems with no symmetry 
for d = 1 (cf. Table Q]). We proceed with a more system- 
atic approach. 



CLASSIFICATION PRINCIPLES 

Concrete mathematical problems may be formulated for 
Dirac operators, band insulators, or more general sys- 
tems. Let us set up the framework. We need to define a 
set of admissible Hamiltonians and some equivalence re- 
lation between them; the equivalence classes may then be 
called "phases". Continuous deformation, or homotopy is 
part of the equivalence definition, but it is not sufficient 
for a nice classification. A key idea in /^-theory is that 
of stable equivalence: when comparing two objects, X' 
and X", it is allowed to augment them by some object Y. 
We generally augment by a trivial system, i.e., a set of 
local, disjoint modes, like inner atomic shells. This cor- 
responds to adding an extra fiat band on an insulator. It 
may be the case that two systems cannot be continuously 
deformed one to the other, but such a deformation be- 
comes possible after the augmentation. Thus, the topo- 
logical classification of band insulators with an unlim- 
ited number of bands is simpler than in the case of two 
bands! Likewise, it is easier to classify Dirac operators if 
we do not impose any restriction on the size of gamma- 
matrices. The final twist is that /^-theory deals with dif- 
ferences between objects rather than objects themselves. 
Thus, we consider one phase relative to another. 

We now give exact definitions for d = (meaning 
that the system is viewed as a single blob). The simplest 
case is where the particle number is conserved, but there 
are no other symmetries. A general free-fermion has this 
form: 

H = ^X jk a)a k: (8) 

j-k 

where X — (Xj^) is some Hermitian matrix representing 
electron hopping. Since we are interested in gapped sys- 
tems, let us require that the eigenvalues of X are bounded 
from both sides, e.g., A < |e 7 | < E max . The following 
condition is slightly more convenient: 

a<e]<a~ l , (9) 

where a < 1 is some constant. This class of matrices is 
denoted by Co(oc), and the corresponding Hamiltonians 
are called admissible. (Some locality condition will be 
needed in higher dimensions, but for d = 0, this is it.) 

The "spectral flattening" transformation, X i— > X = 
sgnX reduces admissible matrices to a special form, 
where all positive eigenvalues are replaced by +1, all 
negative eigenvalues are replaced by —1, and the eigen- 
vectors are preserved. (The matrix element X^ is, essen- 
tially, the equal-time Green function.) Such special ma- 
trices constitute the set 

P)(l)= U U(n)/(U(k)xU(n-k)), (10) 

0<k<n 

where n and k are the matrix size and the numbers of — 1 
eigenvalues, respectively. 



We write X' w X" (orX" w X" to be precise) if X' and 
X" are homotopic, i.e., can be connected by a continuous 
path within the matrix set Cq((x). It is easy to see that 
two matrices are homotopic if and only if they agree in 
size and have the same number of negative eigenvalues. 
For families of matrices, i.e., continuous functions from 
some parameter space A to Co(a), the homotopy clas- 
sification is more interesting. For example, consider an 
integer quantum Hall system on a torus. The boundary 
conditions are described by two phases (mod27r), there- 
fore the parameter space is also a torus. This family of 
Hamiltonians is characterized by a nontrivial invariant, 
the first Chern number II 2 811 . 

It is clear that Cq(oc) can be contracted within itself to 
Co ( 1 ) since we can interpolate between the identity map 
and the spectral flattening: X i— > ft{X), where t £ [0, 1], 
/o(x) = x, /] (x) — sgnx, and the function f is applied to 
the eigenvalues of Hermitian matrix X without changing 
the eigenvalues. Thus, Co(tt) is homotopy equivalent to 
Co(l), and we may use the latter set for the purpose of 
topological classification. 

Let us consider this example (where X is a single 
matrix or a continuous function of some parameters): 

The actual homotopy is Y, = cos(tn/2)Yo + sm(tn/2)Y\. 
Note that Y, 2 = 1 since T 2 = if = 1 and Y Q Y { = -Y { Y . 
Furthermore, Y\ is homotopic to the matrix that consists 
of o z blocks on the diagonal; such matrices will be re- 
garded as trivial. This example shows that any admissi- 
ble system (X) is effectively canceled by its particle-hole 
conjugate (—X), resulting in a trivial system. That is al- 
ways true for free-fermion Hamiltonians, with any sym- 
metry, in any dimension. 

Equivalence between admissible matrices is defined as 
follows: 

X' ~ X" if X' © Y « X" © Y for some Y, (12) 

where © means building a larger matrix from two di- 
agonal blocks. Without loss of generality, we may as- 
sume that Y is trivial. Indeed, if X' © Y w X" © Y, then 
X' © Y © (-Y) w X" © Y © (-F), and we have seen that 
Y © (-Y) is homotopic to a trivial matrix. 

The difference class d(A,B) of two same-sized matri- 
ces is represented by the pair (A , B) up to this equivalence 
relation: 

(A',B') ~ (A",B") if A'ffifi"~A"ffifi'. (13) 

Note that the the matrix sizes in different pairs need not 
be the same. Since (A,B) - (A © (-B), B © (-B)), it is 
sufficient to consider pairs where the second matrix is 
trivial. Thus, the equivalence class of (A,B) is given by a 



single integer, k = k(A) — k(B), where k( - ■■) denotes the 
number of negative eigenvalues. Since B is trivial, k(B) 
equals half the matrix size, n = 2s. Hence, k(A) = s + k. 

To characterize the difference between two families 
of matrices parametrized by A, one needs to consider 
functions from A to the classifying space Co: El 

C = I J lira U (2s) /(U (s + k) xU(s-k)). 

It is the same space as in Table [2] The Abelian group 
of difference classes (= homotopy classes of functions 
A -> C ) is denoted by #°( A ) = n( A ,Co). 



SYMMETRIES AND CLIFFORD 
ALGEBRAS 

In this section, we complete the d — classification. 
Since the particle number is not generally conserved, 
we will use the Hamiltonian Ha given by a real skew- 
symmetric matrix A (see Eq. (O). To generalize some 
arguments of the previous section, let us also define the 
trivial Hamiltonian: H tI i v = — j) = Hq, where 



Q 



\ 



\ 



■J 



(14) 



The eigenvalues of A come in pairs (+z£/, — z£y), 
where e ; - are positive and satisfy inequality ©. Replac- 
ing A with A = — /sgn(;A) takes £, to 1. The matrix A can 
be represented as SQS~ l , where S £ 0(2n). However, 
this representation is not unique since S can be multiplied 
on the right by any orthogonal matrix that commutes with 
Q. Such matrices form a subgroup of 0(2n) that may be 
identified with U(n). Thus, the set of matrices A (i.e., real 
skew-symmetric matrices with ±z eigenvalues) is equal 
to 0(2n) /U(n). Let us take the n — > °° limit by identify- 
ing A with A © Q (where the size of Q can be any even 
number). The result is listed in Table|2]as the classifying 
space R2: 

R 2 = lim 0(2n)/U(n). 

n — >°o 

The set R2 has two connected components, which are dis- 
tinguished by the value of sgn(PfA) = PfA = detS = ±1. 
The physical meaning of this invariant is the fermionic 



parity (— \) N in the ground state, where N 
the particle number. Note that A' is conserved (mod 2) 



3 Here lim m ^„o is a so-called direct limit: the unitary cosets for smaller 
m are mapped into ones for larger m. 



The condition that A' is conserved as an integer is 
equivalent to a t/(l) symmetry. In this case, the creation- 
annihilation expression of Hamiltonian (f2]i should not 
contain terms like djd^ or d\d^. This is a good point to 
note that the approach based on free-fermion Hamilto- 
nians is fundamentally incomplete since it cannot distin- 
guish between the full U(l) group and its Z4 subgroup, 
which is generated by the transformation dj 1— > idj. Let 
us assume for a moment that the actual symmetry is Z4. 
Then terms like aifl2fl3fl4 are allowed in principle, but 
not in a free-fermion Hamiltonian. Therefore topologi- 
cal invariants of noninteracting systems may not be pre- 
served in the presence of interactions. In the following 
example, the number of particle-occupied states changes 
by 4 by a continuous path through an interacting phase: 



H(t) = cos(7W 



4 

E 

.7=1 



d\d j 
j J 



- sin(7Tf) (aifl2fl3fl4 + h.c). 



Note that the ground state remains non-degenerate for all 
values of t. On the other hand, a homotopy like that is 
only possible if the interaction term exceeds the energy 
gap at some point. Thus, the noninteracting topological 
classification is generally stable to weak interactions, but 
not to strong ones. In the U(l) case, it is absolutely stable 
though (at least, for d — 0). We now set this discussion 
aside and proceed with the noninteracting case. 

It is easy to see that the Hamiltonian (fJJ is U(l) 
invariant if and only if the matrix A commutes with 
Q (see Eq. (fT4]>). Another possible symmetry is time- 
reversal invariance. It can be expressed by an antiunitary 
operator T acting in the Fock space; this action is defined 
as follows: 
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Converting , . . . into ci,C2,C3,C4, . . ., we 

obtain a relation of the form Tc m T~ l = Y,i Tinfih where 
the matrix T consists of 4 x 4 blocks: 
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(15) 



The f-invariance of the Hamiltonian is equivalent to the 
condition TA = —AT. 

Let us describe a common algebraic structure that 
is applicable to three symmetry types: no symmetry, T 
only, and T and Q. First, note these identities: 



Q 2 



1. 



TQ=-QT. 



(16) 



It is convenient to introduce some new notation: e\=T, 
<?2 = QJ . Note that ei anticommutes with A if both 
the T and Q symmetries are present. Let us also use 
A = — ;sgn(/A) instead of A. Then we have the following 
characterization: 

No symmetry: A 2 = — 1 ; 
Tonly: e 2 =A 2 = -l, e\A = -Ae\; 
T and Q: 



e 2 = ej=A 2 = -1, 



-e2 e l, : 



-Aej (j = 1,2). 



The pattern is pretty obvious. We have p predefined ma- 
trices e\, . . . ,e p (p = 0, 1,2) satisfying Clifford algebra 
relations (see exact definition below) and look for all pos- 
sible choices of another Clifford generator e p+ \ = A. 

The (real) Clifford algebra Ciff p<? is generated by 
elements ei, . . -,e p+q satisfying these relations^ 



s? = = e 2 = -1 



" 2 , 
'P+l 



'P+q 



= 1, 



= -e*«7 U^k). 



(17) 



All Clifford algebras can be described in terms of the 3 
simple algebras with real coefficients: R (real numbers), 
C (complex numbers), and H (quaternions). For exam- 
ple, Ciff ' 1 is isomorphic to R©R since it consists of 
linear combinations of two complementary projectors, 
i(l ±ex). The algebra Ciff 10 can be identified with C 
by mapping the negative generator e\ to i. Furthermore, 
Ciff 0,2 Ciff 1 ' 1 ^ R(2) (the algebra of real 2x2 matri- 
ces where the Clifford generators are mapped to O z , O x 
or to a z , i<J y , respectively), and Ciff 2 = H. For more 
details on Clifford algebras and their use in ^-theory, see 
Refs. HQ. 

In the problem at hand, the Clifford generators act in 
the mode spaced Thus, we deal with Clifford algebra 
representations such that e\,. 



,e p and e p+u . 



'p+q 



are represented by real skew-symmetric and real sym- 
metric matrices, respectively. To classify free-fermion 
Hamiltonians, we consider representations of Ciff p+ 10 
with fixed action of e\, . . . ,e p ; we call that the "Clifford 
extension problem with p negative generators". 

For technical reasons, it is convenient to reformu- 
late the problem in terms of positive generators. To 
this end, we will employ the isomorphism Ciff' p+2 = 
Ciff' ®R(2), which may be defined as follows: 



ej i ► e'j (g> (io y ) for j = 1 , . 
e p+ i^I®G z , e p+ 2i->I( 



(18) 



An alternative notation is also used, where the positive generators 
(e 2 j = 1) are listed first and the parameters p and q are swapped. 
5 In comparison, the Majorana operators c; generate a (complex) Clif- 
ford algebra acting in the Fock space. 



Representations of the algebra stf <S> M(«) (for any g/) 
have very simple structure, namely, E ®R", where the 
first factor is some representation of and the second 
comes with the standard action of the matrix algebra 
M(n). Thus, s>/ and srf ®R(n) have the same representa- 
tion theory (i.e., their representations are in a natural one- 
to-one correspondence); such algebras are called Morita 
equivalent. Up to Morita equivalence, Ciff'^ only de- 
pends on p — q mod 8. 

We conclude that the classification of free-fermion 
Hamiltonians with p negative Clifford symmetries is 
equivalent to the extension problem with q = p + 2 pos- 
itive generators. That is to say, we need to find all possi- 
ble actions of e q +\ (e 2 +l = 1) if the action of e\ , . . . ,e q 
is fixed. In ^-theory, the problem is formulated in terms 
of difference objects (E,F,w), where E, F are represen- 
tations of Ciff <?+1 and w is a linear orthogonal map 
that identifies them as Ciff ' 9 representations, see 1251 . 
Without loss of generality, we may fix F to be a sum 
of several copies of the regular representation (which 
corresponds to a trivial Hamiltonian) and w the iden- 
tity map. Such difference objects form the classifying 
space R q (see Table The Abelian group of equiv- 
alence classes of difference objects parametrized by A 
is A^^A) = 7t(A,R c/ ). It is isomorphic to the conven- 
tional real ^T-group ^^(A), which is also denoted by 
KO~ q (A). In the special case where A = pt (a single 
point), we get K^ q (pt) = no{R c ,)- 



CLASSIFICATION FOR ARBITRARY d 

We begin with a short summary, focusing on the symme- 
try classes that correspond to real ^-theory. It is natural 
to distinguish three cases: 

1. Continuous free-fermion Hamiltonians are classi- 
fied by /T R <? (5 i/ ) = 7iQ{R q _ c [), where S d represents 
the momentum space (see below). Sufficient insight 
can be gained by considering Dirac operators. This 
setting is actually more general than one might ex- 
pect: gapped Hamiltonians in the momentum space 
are topologically equivalent to nondegenerate mass 
terms that anticommute with a fixed Dirac operator. 
Long-range disorder may be described by textures 
of the mass term varying in space, i.e., continuous 
functions M : W 1 — > R q _ d . 

2. Band insulators are characterized by the momen- 
tum space T d , hence the classification is given by 
K~ q (t d ). This Abelian group includes 7to(Rq-d) as 
a direct summand, but there is some extra piece (cf. 
"weak topological insulators"). 

3. Arbitrary local discrete systems under the energy 
gap or localization condition. ("Local" means that 



the electron hopping is short-ranged. The gap con- 
dition is stronger than the localization, but the prob- 
lem for the localized case can be reduced to that for 
the gapped case.) Realizations of short-range disor- 
der fall into this category. The classification of such 
general systems is exactly the same as for Dirac op- 
erators, due to the following 

Theorem: Any gapped local free-fermion Hamilto- 
nian in W l is equivalent to a texture. 
(That is the key technical result, but I cannot ex- 
plain it in any detail in such a short note.) Discrete 
systems on a compact metric space L are classified 
by the /f-homology group Kf(L). 

Continuous systems and Dirac operators 

The Hamiltonian of a translationally invariant systems 
can be written in the momentum representation: 

where j and k refer to particle flavors. The matrix A (p) is 
skew-Hermitian but not real; it rather satisfies the condi- 
tion A ; ^(p)* = Ajj l (—'p). By abuse of terminology, such 
matrix-valued functions are called "functions from Mr 
to real skew-symmetric matrices", where R'' is the usual 
Euclidean space with the involution p <-> p (cf. [29*]). 
The symmetry is defined by some Clifford generators 
represented by real matrices whose action does not de- 
pend on p. As described in the previous section, we can 
turn negative generators to positive and replace A(p) by 
another Clifford generator e^+^p). While the matrices 
<?i, . . . ,e q are real symmetric, e q+ \ is Hermitian and sat- 
isfies the condition e q+ \ (p)* = e q+ \ (— p). 

A reasonable classification can be developed when the 
asymptotics of A(p) for |p| — > °° is fixed. We may iden- 
tify the infinity in the momentum space with the bound- 
ary of a large ball, dB d . Thus, the difference between two 
phases may be characterized by an element of the relative 
/f-group 

K^\B d , dB d ) = K^(S d ) <* 7Co(R q - d ). (20) 
Here we have used the isomorphism ll25ll 

K^(X)=K^(S r X) (r = mod 8), (21) 
and the (1, 1) periodicity [29]: 

K^(SSX)^K^(X), (22) 

where S denotes the suspension. 

The group Tlo{Rq-d) — (pt) on the right-hand side 
of Eq. ( f20l > has a concrete physical interpretation. It 



classifies the nondegenerate mass terms M in the real 
self-adjoint Dirac operator D = Y.a7ad a +M, where y a 
are skew-symmetric, M is symmetric, and 

YaYb + YbJa = S ab , YcM = -My a . (23) 

(Replacing M with M = sgnM, we can achieve that 
M 2 = 1 .) In addition, we assume that Yi > • • • > Yd ar >d M 
anticommute with the symmetry generators e\,...,e q . 
Thus, the gamma-matrices play the role of Clifford sym- 
metries with opposite sign; they effectively cancel the ac- 
tual symmetries. Note that those new "symmetries" do 
not entail any conservation laws. Our argument only im- 
plies that any continuous spectrum is equivalent (up to 
an augmentation and homotopy) to a Dirac spectrum that 
has the additional symmetries. 



Discrete systems 

Let us consider the Hamiltonian ©, where each mode 
j is associated with a site, or point ry in the real space. 
There may be several modes per site; symmetries (if 
any) act independently on each site. We assume that the 
Hamiltonian is r-local (i.e., A,* = if the distance be- 
tween Tj and r# is greater than r) and that it is a-gapped 
(i.e., the eigenvalues Ej of iA satisfy inequality (0). Un- 
der these conditions, the matrix element A decays very 
fast as the distance between j and k goes to infinity, 
which is a sign of localization. Conversely, if we start 
with the matrix A (such that A 2 = — 1 and A » decays 
fast enough) and set all the elements for |r; — r^j > r' to 
zero, we will obtain a gapped local matrix. Both transfor- 
mations can be done continuously, which roughly shows 
that the set of localizing Hamiltonians is contractible 
within itself to the set of gapped Hamiltonians (up to a 
change of controlling parameters). Thus, we may stick 
with the gapped case without any loss of generality. 

Using the standard trick, we replace A with a real sym- 
metric matrix X that is r-local, a-gapped, and anticom- 
mutes with q positive Clifford symmetries. The above- 
mentioned theorem pertains to such matrices. Here, we 
only discuss it at the physical level. The texture corre- 
sponding to the matrix X is constructed algorithmically, 
albeit in a contrived fashion. The procedure is local, with 
a characteristic radius / = cr, where c depends on d and 
a. The number of Dirac modes needed is l d (for local- 
ized systems, it's the localization volume). To calculate 
M(r), we only look at the /-neighborhood of point r, and 
M doesn't vary much at distances smaller than /. We may 
then discretize the Dirac operator on a fine grid, with y a 
and M rescaled properly so as to keep the a parameter 
fixed. Thus, we obtain an r'-local, a-gapped matrix X' , 
where r 1 is arbitrary small. The equivalence between X 
and X' involves an augmentation and a homotopy, where 



r may increase by a constant factor before it shrinks 
down to /. 

This theorem implies that the boundary between two 
phases must carry some gapless modes. Indeed, each 
phase may be characterized by the mass term M(r) com- 
puted at any point away from the boundary. Since the 
phases are different, the two mass terms, M(r\), M{r2) 
belong to different connected components of the classi- 
fying space. But if the boundary between the phases were 
gapped, we could make the whole system into a continu- 
ous texture, and thus M(r\) and M(r2) would belong to 
the same component — a contradiction. 

A gapped local system on a compact metric space 
L (say, a manifold with or without boundary) is char- 
acterized by a /f-homology class t, G K^(L\ where 
q is defined (mod8). ^-homology (see e.g. IB0I1 ) and 
the related noncommutative geometry i3lll are advanced 
subjects, but the basic intuition is rather simple. Let 
us consider systems with no symmetry (q = 2) on the 
two-dimensional torus, L = T 2 . Such systems are triv- 
ially characterized by the number of fermions in the 
ground state, <i;o G Z2. Now imagine a closed Majorana 
chain winding around the torus. It defines a homology 
class ^1 G //i(T 2 ;Z2), which is a topological invariant 
for gapped local systems. It can be measured by cut- 
ting the torus along some cycle c and counting edge 
modes (mod2). Or one can flip the sign of all matrix 
elements A,j spanning across the cut and see how ^0 
changes: (c) = <^o(+) <^o(— )■ If the torus is filled with a 
p x + ipy superconductor, the system has a nontrivial two- 
dimensional invariant, ^2 G Z. But if ^2 is odd, then the 
properties of the ID invariant change: %\ is not a homol- 
ogy class, but rather, a spin structure. Indeed, 

&(++)&)(+-)§>(-+)§>(—) = (24) 
where ± refers to the sign of matrix elements across 
two basis cycles. In general, the definition of low- 
dimensional invariants (except in dimension 0) depends 
on the higher-dimensional ones@The /^-homology class 
includes all. 



Band insulators 

The main difference from continuous systems is that 
the momentum space is T d . Since there is no need to fix 
the spectrum at infinity, the classification is given by the 
absolute ^-group K R ^(T d ) = K R q (T d ). The band struc- 
ture analysis in Refs. Q El 01 an d others offers a con- 
crete view of that group in certain cases. Unfortunately, 



The term "invariant" is used in a sloppy way, but one can rigorously 
define the range of b, s , assuming that i = ... = §j = 0. It is the 
Abelian group E™ q _ s of the ST-homology spectral sequence. 



the momentum space picture is non very intuitive. To un- 
derstand and calculate the group K R g (T d ), we relate it 
to Zf-homology of the real-space torus by means of the 
Baum-Connes isomorphism for II 1 (a -theory analogue 
of the Fourier transform). Then we apply the Poincare 
duality. Thus, 



K R \f d )=Kf(T d )=Ki-«(T d ) 
= K (R q - d )®K d R q (T d ) 



(25) 



The first term is the same as before, but the last one is 
new. It further splits, though not canonically: 



d-l 



(26) 



For 3D T-invariant insulators, i.e., d = 3, q = 4, we get: 



K R \T)^Z®3Z 2 . 



(27) 



The Z term is the number of (Kramers degenerate) va- 
lence bands, whereas 3Z2 pertains to "weak topological 
insulators". 



THE EFFECT OF INTERACTION 

Topological properties of gapped local free-fermion sys- 
tems are mostly understood. The big open question is 
how the classification is changed by interactions, e.g., 
whether different free-fermion phases can be deformed 
one to another through an interacting phase without clos- 
ing the gap. In some cases, e.g., the integer quantum Hall 
effect and chiral 2D superconductors, the topological in- 
variants are related to physical properties that are well- 
defined in the presence of interactions (namely, the Hall 
conductivity and the chiral central charge, which deter- 
mines the edge energy current 11321 Il4ll ). The Kramers 
degeneracy analysis of vortex-bound states demonstrates 
the stability of 2D topological insulators [2j l33ll and 
(Px+ipy)^ + (p x — iPy)i superconductors 1122 1. 

However, the free-fermion classification is unstable 
for ID systems with the unusual T symmetry: T 2 = 1 
instead of t 2 = (— 1) . For a concrete model, consider 
the Majorana chain and its variations, where T acts on 
odd sites by TcjT~ l — —cj so that terms like icjc^ are 
only allowed between sites of different parity. In the 
free-fermion setting, this symmetry is described by one 
positive Clifford generator, hence p = — 1, q = p + 2 = 
1, and for d = 1 we get a topological invariant k G 
7lo{R q -d) = Z. For example, the usual phase transition 
in 8 parallel Majorana chains is characterized by k = 8. 
But in this particular case, the two phases are actually 
connected through an interacting phase 13411 . 
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